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Abstract] The classificatory schema in some groups of Megachilidae reorganized by means of shape differences in 


terms of the inter- and intra-specific and inter- and intra-generic and subgeneric deformations. The relative warps 
UPGMA trees and PCA graphs obtained from 30 landmarks were shown. Both Anthidium cingulatum Latreille and A . 
florentinun] Fabricius[] should better classify in the same subgenus Anthidium s. str. In the other case under the 


results of shape variation our results supported the hypothesis that Paraanthidium as a subgenus of Anthidium 


Fabricius. On the other hand Icteranthidium Michener and Paraanthidiellum Michener showed enough dissimilarity to 
be kept as separate genera. Our analysis also supported the monophyletic origin of Megachilini which was clustered as 


a separate group by all methods used. 


Key words[] Megachilidae[] Anthidiini[] systematics[] geometric morphometrics] UPGMAL[] generalized procrustes 


analysis[] relative warps 


1 INTRODUCTION 


Resembling most of the other groups in 
entomology{] there is still no consensus among 
researchers about the classification of some groups in 
Megachilidae systematic{] Warncke[] 1980[] Özbek and 
Zanden] 1993[] Gogola[] 1999[] Michener{] 2000[J. 
Possible reasons for this confusion may be the lack of 
sufficient morphological characters or the existence of 
artificial intermediate taxa. Howeverl[] these problems 
were partially solved by means of changes made in the 
basic characters leading to changes in the place of some 
generd] Michener{] 20001. 

The tribe Anthidiini also experiences similar 
systematical confusion mentioned above. The tribe can 
be differed from the others in terms of the following 
features[] Shortness of stigma and pterostigma[] the 
existence of yellow[] white or red maculations on 
integument{] having cleft claws in females or having 
basal articulation[] the existence of the second recurrent 
vein in the distal of the second submarginal cell[ the 
lack of dorsal lamella of metapleura[] the lack of 
preapical transverse flange or carina in T6 of the males 
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[] Banaszak and Romasenko[] 1998[] Michener[] 2000[]. 
The earliest study on the phylogeny was carried out by 
Miillef] 1996[]. Despite the whole variation in the male 
genitaliaL] Warnck¢] 1980[] classified all the genera of 
Anthidiini which are not parasitic in the genus 
Anthidium Fabricius] 1804. However[] the taxa which 
were considered as subgenera in this study were treated 
as genera level by Tkalcu[] 1966[[] Banaszak and 
Romasenkd] 1998[]] Gogold] 19910 1999]. Özbek and 
Zander{] 1993[] regarded Paraanthidium Friese[] 1898 
as a subgenus of Anthidium Fabricius[] 1804 and 
Paranthidiellum Michener[] 1948 as a subgenus of 
Pseudoanthidium Friese[] 1898. Michener and Griswold 

O 1994[] and Michener{] 2000[] viewed Paraanthidium 
Friese] 1898 and Archianthidium Mavromustakis{] 1939 
as subgenera belonging to the genus Trachusa Panzer[] 
1804 and he regarded the Paranthidiellum Michener[] 
1948 as the synonymous of the Pseudoanthidium 
Friese[] 1898. 

The lack of consensus in the systematical studies 
directed several authors toward the seach for alternative 
approaches[] Mayr and Ashlock[] 1991[]. The DNA 
based techniques and cladistic methodology in time 
partially solved some of these problems. The use of 
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geometric morphometrics in understanding the process 
underlying the biodiversity is a relatively new 
methodology. After the big revolutionary change in the 
1990’s[] this approach began to represent an important 
new paradigm for the statistical study of“ shape” in 
biology and other fields of science[] Rohlf] 1999a[] 
Rohlf] 2002[]. In the last ten years[] geometric 
morphometric citations per year have increased 
considerably] Lynch 2004]. Especially the landmark- 
based morphometrics has gained significant support 
among scientists including anatomist¢] Pavlinov[] 20010] 
Hennessy and Stringer[] 2002[] Lockwood et al . [20020 
and entomologist{] Alibert et al .[] 2001] Kligenberg[] 
20030. Unlike the analytical approaches[] the geometric 
one is aimed at comparison of the shapes themselves 
[] Pavlinov[] 20010. Its high level of explanation power 
on the shape and size variation and co-variation among 
the organisms[] provided an alternative solution 
especially for problematic groups in systematics and 
ecology] Aytekin et al .[] 20030. 

The principles of these methods are based on 
capturing the two- or three-dimensional Cartesian 
coordinates of landmarks which are basically the 
homologous points among the individuals that have been 
previously assigned the same names [] Bookstein[] 
1991[]. Differences among individual configurations of 
landmarks can be translated to several mathematical 
functions which fit the differences[] Alibert et al. [O 
2001[]. Then the Procrustes distanc] the square root of 
the sum of the summed squared distances between 
homologous landmarks[] can be used as a metric for 
comparing shapes[] Rohlf] 1999a[]. The registered 
landmark configurations [] k landmarks and m 
dimensions[] can be represented as points in a shape 
space which is of lower dimensionality than the figure 
space since location[] rotation and scale differences have 
been removed] O’Higgins and Jones[] 1999[]. For two- 
dimensional data[] the space of Procrustes registered 
specimens is therefore of dimensionality km-4 

O Bookstein[] 1991[] O’Higgins and Jones|] 1999[] Rohlf{] 
1999a[]. This space is known as Kendall’s shape space 
O Kendall[] 1984] which is however non-Euclidian. But 
it possibly can allow the study of well known 
multivariate analysis[] Sneath and Sokal] 1973[] like 
obtaining a principal components analysis] PCA[ or 
unweighted pair-group method using arithmetic averages 
O UPGMA[] phenogram if it is superimposed orthogonally 
[ or stereographically[] in the tangent plane to Kendall’s 
shape spacd] Dryden and Mardia{] 1998[]. In addition[] 
as these methods preserve the geometrical relationships 
among landmarks[] they also provide a powerful vision 
for visualizing deformations in the shape of the original 
specimen in means of displacement vectors or 
deformation grids as warps[] Rohlf] 1993[] Alibert et 
al .{] 2001[]. Some structures like the insect wings with 


the least degrees of freedom are one of the most 
appropriate structures for such studies [] Pavlinov[] 
2001[]. Together with warps PCA of Procrustes 
residuals and UPGMA phenograms can provide an 
excellent combination of techniques for the two purposes 
of such studies as we made for Anthidiini here[] First to 
detect and then to describe the differences among taxa 


[] Lockwood et al .[] 2002[]. 
2 MATERIALS AND METHODS 


2.1 Study area and sampling 
92 specimens were used from the Megachilidae 
collection in Hacettepe University which were sampled 
between 1999 and 2003 from different eco-regions of 
Middle Anatolia{] Fig. 1[]. The study areal] totally 
178.080 km’ [] has generally a steppe character with 
Astragalus sp.[] Acantholimon sp.[] Centaurea sp. 
Echium sp.[] and Anchusa sp. There are mountains 
and hill$] altitude ca. 1 800 ml] that covered mainly by 
Pinus nigra Arn. ssp. pallasiana{] Lamb[] Holmboe[] 
different species of Quercus sp.[] Picea sp. and 
Juniperus sp. with Cedrus libani A. Richard] Davis[] 
1968 — 1985[]. Agriculture and urbanization with high 
industry is also common in most places where 
collections were made. Right front wings[] right hind 
wings and headd] from frontal view[] of the females of 
24 Anthidium cingulatum Latreille 18090 12 
Anthidium florentinum [|] Fabricius] 1775 O 15 
Icteranthidium laterale Latreille[] 18090] 15 Megachile 
maritima [] Kirby 1802 [1] 16 Paraanthidium 
interruptum [] Fabricius[] 1781 [ and 10 
Paraanthidiellum lituratum [| Panzer] 1801 [] were 
prepared for morphometric analysis. Megachile 
maritima which is from Megachilini was used as control- 
group for the analysid] not as reference[]. Megachilini 
were preferred because the tribe Osmiini was understood 
to be paraphyletic with respect to Megachilini] Engel[ 
1999[]. Diagnose of the species and nomenclature 
followed that proposed by Banaszak and Ramasenko 
O 1998] and Warnckd] 1980[]. All the specimens were 
screened at the time they were collected for the presence 
of known ecto- and/or endo-parasites to shield the 
morphometric data from possible traumatic variations 
[O Mayr and Ashlock[] 1991[] Aytekin et al .{] 20020. 
The body parts were then removed from each specimen 
with forceps and mounted in entellane on labelled 
slides. All slides were photographed using a Leica MZ- 
7.4] r| stereoscopic zoom dissection microscope with a 
DC-300 digital camera system[] digitized[] and archived 
with special information codes added. From the 92 
specimens some of them were eliminated because of the 
problems during preparation procedure and finally 92 
specimens for front wings[] 91[] one specimen from P. 
interruptum is ignored[] specimens for hind wings and 
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specimens were digitized twice [] Arnqvist and 
Martensson[] 1998 []. The second session of 
measurement was conducted after the specimens had 
been removed and replaced under the microscope in 
order to take the positioning error into account 
O Arnqvist and Martensson[] 1998[] Alibert et al .U 
2001[]. No analogous systems were used during the 
whole procedure. 


3 RESULTS 


3.1 Relative warps 
By using an orthogonal alignment projection 
method data obtained from front wings[] hind wings and 
heads were analyzed by means of the first two relative 
warps. For the front wings[] relative positions of the 
average configurations of the species A. cingulatum{] 
A. florentinum and P. interruptum were clustered 
together where P. lituratum and I. laterale formed 
another group while M. maritima separated from rest of 
all in the shape space defined by the first two relative 
warps] x = 19 y = 20 a = ODO Fig. 50. The 
deformation grids of the front wing landmark 
configuration showed that the out-group M. maritima 
have a typical shape difference than the other taxa 
O Fig. 60.0 Eigenvalue was calculated as 4.788 for the 
first relative warp[]. Whereas Anthidium spp. and 
Paraanthidium interruptum have similar front wing 
shapes when they are superimposed on the mean shape 
configuration. Same methodology also used for the hind 
wings and heads. We have got similar configuration for 
Anthidium spp. and Paraanthidium interruptum on the 
hinwings. Relative positions of the 
configurations of the species A. 


average 
cingulatum[] A. 
florentinum and P. interruptum were lineerly arranged 


where P. lituratum and M. maritima formed another 
group while J. laterale separated from rest of all in the 
shape space defined by the first two relative warp$] x = 
10y = 20 a = OW Fig. 70. The deformation of the 
hindwings for A. florentinum and P. interruptum were 
more similar especially on the landmarks 5 and 4] Fig. 
$C Eigenvalue was calculated as 2.107 for the first 
relative warp[]. This is also determined on the head 
shape. Relative positions of the average configurations 
of the species A. cingulatum] A. florentinum{] I. 
laterale and P. interruptum were clustered together 
where P. lituratum and M. maritima separated from 
rest of all in the shape space defined by the first two 
relative warp{] x = 10 y = 20) a = OL] Fig. 90. The 
similar orientation of the landmarks 10 3[] 4[] 5 and 6 
showed typically for Anthidium spp. and 
Paraanthidiunt] Fig. 100.[] Eigenvalue was calculated 
as 2.330 for the first relative warpL]. 
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Fig. 5- Relative positions of the average configurations of the 
species 1] A. cingulatum{] 2) A. florentinum{] 30 I. 
laterale] 4] M. maritimal] 5U P. interruptum{] 6U P. 
lituratum[|] for Anthidiini and Megachilini in the shape space 
defined by the first two relative warp$] x = 10 y = 20 a =O. 


Data from front wings. Filled points indicate the species. 
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Fig. 6 Deformation grids of the six species of Anthidiini and Megachilini used in the present study superimposed on the reference 
configuration] x = 1 y =2[] a = OL] in relative warps. Data from front wings. Open circles with numbers indicate the landmarks. 
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Fig. 7 Relative positions of the average configurations of the 
species] 10 A. cingulatum] 20 A. florentinumi 30 1. 
laterale] 4|] M. maritima|] 5U P. interruptum|] 6] P. 
lituratum[|] for Anthidiini and Megachilini in the shape space 
defined by the first two relative warp$] x = 10 y = 20 a = OÙ. 
Data from hind wings. Filled points indicate the species. 
































3.2 PCA 
When a PCA is conducted by using the data 
obtained from the shape of the front wings[] hind wings 


and heads we reached similar results. The species were 
ordered along the first two principal components similar 
to front wings[] Fig. 11[] and hind wings{] Fig. 120 
data. The PCA of head data gave relatively sparse 
distribution than the ones obtained from the wing$] Fig. 
130. 
3.3 UPGMA 
The UPGMA phenogram formed by SAHN 
clustering obtained from the data of whole landmarks 
O 60 characters[] used in the study showed distinct group 
patterns for all taxa[] Sneath and Sokal] 19730. 80 
specimens were randomly used for the cluster formation . 
The control-group M. maritima clustered separately 
from the rest of the taxa. I. laterale and P. lituratum 
also clustered as a distinct group which indicates not 
only the species difference but also enough for 
identification of the genus level. A. florentinum[] A. 
cingulatum and P. interruptum showed enough distinct 
patterns for being different species group but showed a 
mixed structure for a different genus category level 


O Fig. 140. 





M. maritima 
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P. interruptum 





Deformation grids of the six species of Anthidiini and Megachilini used in the present study superimposed on the reference 


configuration] x = 1[ y = 2[] a = OL] in relative warps. Data from hind wings. Open circles with numbers indicate the landmarks. 
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Fig. 9 Relative positions of the average configurations of the 
species 1] A. cingulatum|] 2) A. florentinum|] 30 I. 
laterale] 4] M. maritima|] 5U P. interruptum|] 6Ų P. 
lituratum[|] for Anthidiini and Megachilini in the shape space 
defined by the first two relative warp$] x = 10 y = 20 a =O). 
Data from heads. Filled points indicate the species. 






































4 DISCUSSION 


The research field of geometric morphometrics is a 
diversity of morphological structures which may be 
uncertain individual variation|] differences among 
discrete groups like taxa[] sexes[] ages[] insect castes[] 
biomorphs] ecotones and others[] Pavlinov[] 2001[] or 
sometimes fluctuating asymmetry] Klingenberg[] 20031. 
Therefore we hypothesized that the classificatory schema 
in some groups of Megachilidae would be reorganized by 
means of shape differences in terms of the inter and 
intra-specific and inter and intra-generic and subgeneric 
deformations . 
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Fig. 10 Deformation grids of the six species of Anthidiini and Megachilini used in the present study superimposed 


2[] a = OL] in relative warps. Data from heads. 


Open circles with numbers indicate the landmarks. 
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Distribution of the six different species of Anthidiini and Megachilini. Principal component analysis 


Fig. 11 


conducted for the 16 landmarks digitized from the front wings. 
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Fig. 12 Distribution of the six different species of Anthidiini and Megachilini. Principal component analysis conducted 
for the 6 landmarks digitized from the hind wings. 
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Fig. 13 Distribution of the six different species of Anthidiini and Megachilini. Principal component analysis 
conducted for the 8 landmarks digitized from the heads. 
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The results obtained from the relative warps of a 
total of 30 landmarks from three main body parts showed 
that A. cingulatum|] A. florentinum and P. 
interruptum did not show enough differences in the 
bending energies of the grid deformations. Therefore it 
is better not to classify those in separate genera. The 
inter-specific variations were basically in landmarks 10 
and 11 in front wings. This information can indicate a 
close irrelativeness but clearly need additional data 
especially from DNA sequence. The general inter- 
specific variations were in landmarks 2[]3 and 5 in hind 
wings while and 4[] 5 and 6 in head shape. 

The patterns of covariance in the landmark 
coordinates after Procrustes superimposition were 
studied with PCA [ Dryden and Mardia{] 19980 
Klingenberg[] 2003[]. The extracted first two PC’s from 
our analysis had a clear geometric explanation in 
Anthidiini. Especially the data obtained from the wings 
supported our hypothesis. Relatively sparse distribution 
pattern from the head data would possibly due to the 
landmark types[] measurement error and most degrees of 
freedon{] Readers unfamiliar with the techniques may 
be directed to Bookstein[] 1991[] and Zelditch et al. 

O 2004[] for details[]. The UPGMA phenogram also 
clearly gave the same schema. 

These results are however different from those of 
proposed by Warncke] 1980[]. Both A. cingulatum 
and A. florentinum showed separate species character 
but should better classify in same subgenus as 
Anthidium s. str. In the other case P. interruptum 
should better not classified as a separate genus. Under 
the results of shape variation we would propose 
Paraanthidium as a subgenus of Anthidium Fabricius[] 
1804. On the other hand [Icteranthidium and 
Paraanthidiellum showed enough dissimilarity to be 
kept as separate genera. This classification supported 
the approach of Özbek and Zanden[{] 1993{]. Our 
analysis also proved the monophyletic origin of 
Megachilini which was clustered as a separate group by 


all methods used. 
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